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Abstract. This paper is concerned with the local bifurcation analysis 
around typical singularities of piecewise smooth planar dynamical sys- 
tems. Three— parameter families of a class of non— smooth vector fields 
are studied and the tridimensional bifurcation diagrams are exhibited. 
Our main results describe the unfolding of the so called fold — cusp 
singularity by means of the variation of 3 parameters. 



1. Introduction 

NSDS's have become certainly one of the common frontiers between Ma- 
thematics and Physics or Engineering. Problems involving impact or fric- 
tion are piecewise— smooth, as are many control systems with thresholds. 
Many authors have contributed to the study of Filippov systems (see for 
instance [7] and |10j). One of the starting points for a systematic approach 
to the geometric and qualitative analysis of non— smooth dynamical systems 
(NSDS's, for short) is |13| . on smooth systems in 2— dimensional manifolds 
with boundary. The generic singularities that appear in NSDS's, to the best 
of our knowledge, were first studied in [15]. Bifurcations and related pro- 
blems involving or not sliding regions were studied in papers like [6l [HJ [H [3] . 
The classification of codimension— 1 local and some global bifurcations for 
planar systems was given in [11] . In [9] codimension— 2 singularities were 
discussed and it was shown how to construct the homeomorphisms which 
lead to topological equivalences between two NSDS's when the discontinuity 
set is a planar smooth curve. See [16] or [1] for a survey on NSDS's and 
references there in. 

The specific topic addressed in this paper is the qualitative analysis of 
fold— cusp singularities of NSDS's, where a fold and a cusp coincide. More- 
over, the bifurcation diagrams are exhibited. 

Specifically, we distinguish the following cases (see Figure H]): 
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• Unfolding of an invisible fold— cusp singularity: 



(1) 



1 

-x + A 



-x 2 + P-!ji(x,/3,»)) 



if y > 0, 
if y < 0, 



where (A, /3) G (— Ao,Ao) x (— /3o,/3o), with Ao > and /3o > sufficiently 
small and B is a C 2 —bump function such that B(x, /3, fj.) = if /3 < and 

), if x < -y 7 ? or x > 4V/3; 



(2) 



B(x,/3,n) 



if /3 > 0, where 



Bfafl + ffafi), if-^<z<\/A 
B 2 (x,(3) + f(P,v), if^<x<4^p. 

x 2 (208 + 3,S)- 



£ 2 (x,/3) 



128^ V 4^^(176 + 15/3) + /3(688 + 93/3) 
-1 ( (x-4V£) 3 ((x 2 + /3)(-16 + 9/3)- 



48/3 



2x^(16 + 15/3)) 



and 



8/3(128 + 3/3)/iV?(256 + 63/3)/x - (-64 + 45/3)/i 2 - 
^-V2(80 + 3/3)^ + /3~ 1 (— 16 + 9/3)^ 4 ) 

Unfolding of a visible fold— cusp singularity: 

1 

x — A 
1 

-x 2 + /3 

where (A, (3) G (— Ao,Ao) x (— /3o,/3o), with Ao > and /3o > sufficiently 
small. 



(3) 



Y R -- 



if y > o, 
if y < o, 





FIGURE 1. Fold— cusp singularities. Following the notation of [9] , 
a point po S S is a fold— cusp singularity of Z = (X, Y) if it is a 
S— fold point of X and a E— cusp point of Y (see precise definitions 
in Section [2]). 
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1.1. Setting the problem. Denote both, Zr\ s u, in (JTJ) and Zp s \ in ([3]), by 
Z = (X, Y). In short our goal is to study the local dynamics of Z consisting 
of two smooth vector fields X and Y in M? such that on one side of a smooth 
surface S = {y = 0} we take Z = X and on the other side Z = Y. 

In [9] the analysis of the bifurcation diagram of the 2— parameter family 



of NSDS's presenting an invisible fold— cusp singularity is performed. A 
challenging problem is to extend the analysis of [9] in answering the following 
question: Can we find families of NSDS's presenting fold— cusp singularities 
whose dynamics is richer than the family exhibited in [9] ? In this paper 
such an extension has been carried out. By means of the positive answer to 
the previous question, we are able to say that two parameters is not enough 
to explain the birth of some new topological types around Zo,o,o- 

In fact, ours results cover the study done in [9] and we can obtain the 
bifurcation diagram presented in [9] assuming (3 = fi 2 and fi < in ([T]). For 
example, the configuration in Figure [2] is not observed in [9] and is present 
at the bifurcation diagram of ([!]). 





Figure 

2 . Configuration 
nearby Zo,o,o n0 ^ 
observed in [9|. 



FIGURE 3. The local and the 
global bifurcation observed in 
(HJ when /3 > 0, A = y 7 ? and 
[1 = 0. 



We mention two particular situations illustrated in Figure [3] that occur 
in dl]) when (5 > 0. In this resonant configuration we note, simultaneously, 
a two— fold singularity (which is a local phenomenon) and a loop passing 
through the visible S— fold of Y (which is a global phenomenon). 

1.2. Statement of the Main Results. Theorems 1, 2 and 3 pave the way 
for the proof of Theorem A. Theorem B is self contained. 

Theorem 1. If fx = in Equation ([T|) then its bifurcation diagram in the 
(A, j3)— plane contains essentially 17 distinct phase portraits (see Fiaure[23\). 
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It is easy to see that the cases covered by Theorem 1 do not represent 
the full unfolding of the invisible fold— cusp singularity. Because of this, the 
next two theorems are necessary. 

Theorem 2. If < fi < /io in Equation (pQ) then its bifurcation diagram in 
the (A, /?)— plane contains essentially 19 distinct phase portraits (see Figure 

B3p. 

Theorem 3. If — /io < fJ> < in Equation ([I]) then its bifurcation diagram in 
the (A, /?) — plane contains essentially 19 distinct phase portraits (see Figure 

ESP. 

Finally, we are in position to state the main results of the paper. 

Theorem A. The bifurcation diagram of Equation (pQ) exhibits 55 distinct 
cases representing 23 distinct phase portraits (see Figure[2l\). 

Theorem B. The bifurcation diagram of Equation (J3j) exhibits 11 distinct 
phase portraits (see Figure\3^ . 

The paper is organized as follows. In Section [2] we present some basic 
elements on the theory of NSDS's. In Sections (3J H] and [5] we pave the way 
for the proofs of the main results of the paper (Theorems A and B). Section 
[6] is devoted to prove Theorem A and exhibit the Bifurcation Diagram of 
(p} . In Section the proof of Theorem B and the Bifurcation Diagram of 
(J3]) are presented and in Section [8] some concluding remarks are discussed. 
In our paper we follow basically the terminology and the approach of [11] 
or [9] and no one sophisticated tool is needed. 



2. Preliminaries 

Let K C IR 2 be a compact set such that dK is a smooth 1— manifold 
and E C X given by E = / -1 (0), where / : K — > M. is a smooth func- 
tion having G E as a regular value (i.e. Vf(p) / 0, for any p £ /~ 1 (0)) 
such that dK n £ = or dK rti E. Clearly the switching manifold E 
is the separating boundary of the regions E + = {q G K\f(q) > 0} and 
E- = {q € < 0}. 

Designate by % the space of C 1 — vector fields on K endowed with the 
C 1 -topology. Call $7 = Q>(K, f) the space of vector fields Z : K — > R 2 such 
that 



(4) Z(x,y) 



X(x,y), for (x,y)eE+, 
Y(x,y), for (x,y) G E _ , 
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where X = (fi,gi), Y = (/2,<72) are in %. We write Z = (X, Y), which we 
will accept to be multivalued in points of E. We endow with the product 
C 1 — topology. The trajectories of Z are solutions of q = Z(q), which has, 
in general, discontinuous righthand side. The basic results of differential 
equations, in this context, were stated by Filippov in [7j. 

Definition 1. A k— parameter family of elements in Q is a C 1 — mapping, 
with r > 1, 

(: S k — ► 

£» = (£1, £>2, ■ • • , Qk) ^ ^ 
w/iere S* = [— e 1( ei] x [— e 2 , e 2 ] x . . . [— e*, Cfc] £j > 0, i = 1, 2, . . . , k, 
sufficiently small. 

Definition 2. We say that W,W G x defined in open sets U and U , re- 
spectively, are C° — orbitally equivalent if there exists an orientation pre- 
serving homeomorphism h : U — )■ U that sends orbits of W to orbits of W . 
Here, orbit ofW means the image of a solution of x = W(x). 

Definition 3. Two non— smooth vector fields Z = (X,Y), Z = (X,Y) E 
Q(K,f) defined in open sets U, U C K and with switching manifold E are 
E— equivalent if there exists an orientation preserving homeomorphism h : 
U — > U that sends E in E, the orbits of X restrict to U n E + in the orbits 
of X restrict to U (1 E + ; and the orbits ofY restrict to U n E~ in the orbits 
of Y restrict to U n S~ . 

Consider the notation 

X.f{p) = (Vf(p),X(p)) and X\f(p) = (VX i ~ l .f(p),X(p)),i > 2 
where (., .) is the usual inner product in M 2 . 

Remark 1. The vertical dotted lines present in almost all figures of this 
paper represent the points p G K Cl 2 where X.f(p) = or Y.f(p) = 

We distinguish the following regions on the discontinuity set E: 

(i) E c C E is the sewing region if (X.f)(Y.f) > on E c . 

(ii) E 6 C E is the escaping region if (X.f) > and (Y.f) < on E e . 
(hi) E s C E is the sliding region if (X/) < and (Y.f) > on E s . 
Consider Z G f2. The sliding vector field associated to Z is the vector field 

Z s tangent to E s and defined at q G E s by Z s (q) = m — q with m being 
the point of the segment joining q + X(q) and q + Y(q) such that m — q is 
tangent to E s (see Figure HJ. It is clear that if q G E s then q G E e for — Z 
and then we can define the escaping vector field on E e associated to Z by 
Z e = —(—Z) s . In what follows we use the notation Z s for both cases. 

We say that q G E is a E— regular point if 
(i) (X.f(q))(Y.f(q)) > or 
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FIGURE 4. Filippov's convention. 



(ii) (X.f(q)){Y.f(q)) < and Z s (g) ^ (that is q G S e U S s and it is 
not an equilibrium point of 

The points of £ which are not £— regular are called T,— singular. We dis- 
tinguish two subsets in the set of X— singular points: X* and S p . Any q G S p 
is called a pseudo equilibrium of Z and it is characterized by Z^(q) = 0. Any 
g G E* is called a tangential singularity and is characterized by Z^(q) ^ 
and (X.f(q))(Y.f(q)) = (g is a contact point). 

We say that a point pq G X is a J]— fold point of X if X.f(po) = but 
X 2 .f{po) 7^ 0. Moreover, po S S is a visible (respectively invisible) X— fold 
point of X if X.fipo) = and X 2 .f(p ) > (respectively X 2 .f(p ) < 0). 
We say that a point qo G £ is a S— cwsp point of Y if Y./(go) = ^ 2 -/(<7o) — 
and y 3 ./((7o) 7^ 0. Moreover, a S— cusp point qo of F is of kind 1 (respec- 
tively kind 2) if Y 3 .f(qo) > (respectively Y 3 .f(qo) < 0). In particular, 
S— fold and S— cusp points are tangential singularities. 

A pseudo equilibrium q G TP is a S— saddle provided that one of the 
following conditions is satisfied: (i) q G S e and g is an attractor for Z s or 
(ii) 5 6 E s and q is a repeller for Z s . A pseudo equilibrium q G S p is a 
T,— repeller (resp. S— attractor) provided g G S e (respectively g G T, s ) and 
g is a repeller (respectively, attractor) equilibrium point for Z s . 

Given a point g G S c , we denote by r{q) the straight line through q + X(q) 
and g + Y(q). 

Definition 4. T/ie T,— regular points q G S c stzc/i £/ia£ either {X(q),Y(q)} 
is a linearly dependent set or r{q) n X = are called virtual pseudo equi- 
libria. 
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Let us consider a smooth autonomous vector field W defined in an open 
set U. Then we denote its flow by <j>w(t,p). In this way, 

( d 

—(/> w (t,p) = W((p w (t,p)), 
i at 

k (f>w(0,p) = P, 

where t £ I = I(p, W) C R, an interval depending on p G U and W. 
The following definition was stated in [9], pg 1971. 

Definition 5. The local trajectory of a NSDS given by Q is defined as 
follows: 

• For p G S + andp G £~ the trajectory is given by (f>z(t,p) = (j>x(t,p) 
and <pz(t,p) = (f>y(t,p) respectively, where t £ I. 

• For p G S c such that X.f(p) > 0, Y.f(p) > and taking the origin 
of time at p, the trajectory is defined as (pz(t,p) = <pY(t,p) for t G 
Id {t < 0} and 4>z(t,p) = 4>x(t,p) for t G I Ci{t > 0}. For the case 
X-f(p) < and Y.f(p) < the definition is the same reversing time. 

• For p G S e UE s such that Z T '(p) we define <f>z{t,p) = 4>z^{t,p) 
for tel. 

• For p G <9X C U <9E e U dT, s such that the definitions of trajectories 
for points in a full neighborhood of p in S can be extended to p and 
coincide, the trajectory through p is this trajectory. 

• For any other point (pz(t,p) = p for all t G R. This is the case of 
points in <9X C U <9£ e U <9X S which are not regular tangential singula- 
rities and the equilibrium points of X in S +; ofY in £_ and of 

in S s US e . 

Definition 6. The local orbit— arc of the vector field W passing through 
a point p G U is the set Jwip) = {<f>w(t,p) - t E I}. 

Since we are dealing with autonomous systems, from now on we will use 
trajectory and orbit— arc indistinctly when there is no danger of confusion. 

Definition 7. Consider Z = (X, Y) G Q. 

n 

(1) A canard cycle is a closed curve T = ^J<7j composed by the union 

i=i 

of orbit— arcs <jj, i = 1, . . . ,n, o/X| s +, Y\j^- and Z s such that: 

• Either there exists iq C { 1 , . . . , n} with o~i C jx ( respectively 
ai C 7y ) and then there exists j ^ io with aj C 7y U 7^s 
(respectively o~j C jx U 7^e ), or T is composed by a single arc 
o-i ofZ s ; ' 

• the transition between arcs of X and arcs ofY occurs in sewing 
points; 

• the transition between arcs of X (or Y) and arcs of Z^ occurs 
through Ti—fold points or regular points in the escaping or slid- 
ing arc, respecting the orientation. Moreover i/T / S then 
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there exists at least one visible T,—fold point on each connected 
component o/mS. 

(2) A canard cycle T of Z is of: 

• Kind IifT meets £ just in sewing points; 

• Kind IlifT = £; 

• Kind III ifT contains at least one visible H—fold point of Z. 
In Figures and [7] arise canard cycles of kind I, II and III 

respectively. 

(3) A canard cycle T of Z is hyperbolic if one of the following conditions 
are satisfied: 

(i) r is of kind I and r]'(p) ^ 1, where r/ is the first return map 
defined on a segment T with p £ T iti 7; 

(ii) r is of kind II; 

(hi) T is of kind III, WnWnF = and either T n E C S c U S e U S* 
or m E C S C UE S UE*. 






Figure 

5 . Canard 
cycle of 
kind I. 



Figure 

6 . Canard 
cycle of 
kind II. 



Figure 

7. Canard 
cycle of kind 
III. 



Following Theorem 2 of [T7], locally it is possible to consider f(x,y) = y 
and conclude that any X £ E + presenting a E— fold point is C°— orbitally 
equivalent to the normal form Xo(x,y) = (p\,p2x) with p\ = ±1 and 
92 = ±1- 



Following |13j . we can take f(x,y) = y and derive that any Y £ E 
presenting a E— cusp point is C°— orbitally equivalent to the normal form 
Y (x,y) = (p3,Pax 2 ) with ,93 = ±1 and p 4 = ±1. 

Lemma 8. let Z = (X, Y) £ Q presenting a fold— cusp singularity, then Z 
is E— equivalent to the standard form Zq = Zp lt p 2tP3tP4 given by 



(5) 



J Pl;P2,P3,PA 



X 



1", 



P3,p4 



Pi 

P3 
p 4 X 2 



ify>o, 
ify<o 



where pi,P2,P3,P4 = ±1. 
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Observe that the values of pi, i = 1,2,3,4, in Lemma [8] depend on the 
orientation of X and Y. In Subsection 12.31 we prove Lemma El i.e., we ex- 
hibit the homeomorphism that characterizes the equivalence between any 
fold— cusp singularity and the standard form given by ([5]). 

Consider Z™ b ' kl , z^ ls ' k2 E O written in the following standard forms (sim- 
ilar forms were stated in Section 12 of [9]): 



(6) Z { 



yivb 



Vkl 



-x 

-1 

„2 



if y > o, 

if y < 0, and 



(7) 



VVIS 



Vk2 



if y > o, 
if y < 0. 



Note that .X^ 6 presents an invisible S— fold point on its phase portrait, Xq is 
presents a visible X— fold point, Yq 1 presents a X— cusp point of kind 1 and 
Yq 2 presents a X— cusp point of kind 2. Moreover, in ([6]) we made pi = 1 and 
P2 = P3 = Pa = —1 in © and in ([7|) we made p\ = p2 = P3 = 1 and /)4 = — 1 
in ©. For simplicity we restrict our study to the normal forms given above, 
i.e, © and ([7]). All the other choices on the values of pi, i = 1, 2, 3, 4 in (JSJ 
are treated similarly. 

The main problem is to exhibit the bifurcation diagram of Zq' p where 
r = ivb or vis and p = kl or k2. 

In order to detect a larger range of topological behaviors near an invis- 
ible fold— cusp singularity we have to refine the analysis done in [9]. This 
refinement can be obtained adding a bump function on the expression of the 
NSDS. 

Denote 

g 2 (x, /3, p)dx = — - (3x + B(x, (3, p) + c , 

where Co = —2f3yf]3/2> and 52 is the second coordinate of Yg jjU in ([1]). The 
C 1 — bump function B satisfies the following properties when /3 > 0: 

• It has exactly one point of local minimum in the interval (— yT?, 4-^//?) . 
This point is located at xq = yffi. 

• F(3yf]3 + ;u) = (see Figured]). By means of this last property the 
orbit— arc of Yg )M that has a quadratic contact to S at go = (— y 7 ^, 0) 
turns to collide with X at the point q\ = (3\//3 + ^,0). So, the first 
coordinate of q\ is bigger (respectively, smaller) than 3y/]3 as p is 
bigger (respectively, smaller) than 0. 
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3# ,3# V? 






[ 1 1 







FIGURE 8. Graph of B. 



Remark 2. It is worth saying that the parameter fi breaks the strong pro- 
portionality between the roots of g2(x, (3, 0). At the limit value ^q = 0, Z\^^ 
presents distinct topological behaviors for fj, < fj,Q or I 1 > Mo- 
Remark 3. Note that in Equations ([T]) and (|3]) the perturbations considered 
depend only on the variable x. The local geometry of a NSDS presenting a 
cusp— fold singularity becomes rather different if perturbations involving the 
variables x and y are admitted. 

2.1. Global Bifurcation. As said before, the configuration illustrated in 
Figure [3] plays a very important role in our analysis. The configuration of 
this figure is reached from (pQ), by taking (3 > 0, A = yffi and [i = 0. In this 
section we deal with this global phenomenon. 

Emphasizing, let Zq = (Xq,Yq) G Vt having the following properties: 

• The discontinuity set X is represented by f(x, y) = y. 

• Consider Xq = and Yq = (/j, g®) and assume that flip) > 
if p G £+ and / 2 °(p) < if p G E~. 

• go S S is a visible S— fold point of Yq and X.f(qo) > 0. 

• The orbit 7x (<?o) of Xq through qo meets transversally S at a point 
?l- 

• The orbit 7y (gi) of Yb through gi meets tangentially S at go- Call T 
the degenerate canard cycle composed by 7x ( ( 7o) an d 7r (gi)- Let 
M be the compact region in the plane bounded by T. 

2.1.1. Transition Fold Map. As go G X is a visible X— fold point of Yq, 
we may assume (see [17] ) coordinates around go such that the system is 
represented by (x,y) = (— l,x) with g = (0,0). The solutions of this 
differential equation are given by: 

<t>a,b{t) = H + a, ~(t 2 /2) +at + b). 

The orbit— arc 4>q through (0,0) is represented by (poit) = (—t, —t 2 /2). 

Let 5 be a very small positive number. We construct the Transition Map 
£ : L\ — > Lq, from L\ = {(x,y), y = —5, x > V26} to Lq = {(x,0),x > 0}, 
following the orbits of Yq (see Figure [9|). The curve L\ is transverse to 
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Yq at ps = (y/25,—5). Since the solutions 4>& through (x, — 5) G L\ meet 
S = {y = 0} at time i = x ± V^ 2 — 2<5 we obtain that = Va; 2 — 2<5 and 
£ is an homeomorphism. Moreover, ^~ x [x) = \J x 2 + 2(5, i;" 1 is differentiable 
at and (f -1 )'(()) = 0. 




FIGURE 9. First Return Map around the two— fold singularity p. 

2.1.2. First Return Map associated to T. Let gx be the transition map from 
Lq to L2 C £ via Xq— trajectories and gy be the transition map from L2 to 
L\ via Yo— trajectories (see Figure [9]). Observe that the linear part of the 
composition gy o gx is nonzero due to the transversality conditions of the 
problem. For simplicity, let J e = [0, e) x {0} be a small semi-open interval 
of S. 

So the First Return Map of Zq at qo is n{x) = (£0 g Y o gx)(x) for x £ J € . 
Its inverse k _1 is a differentiable map at and satisfies (k~ 1 )'(0) = 0. So, T 
locally repeals the orbits of Zq closed to V and in the interior of V. 

In conclusion, if Z is very close to Zq in Q in such a way that it possesses 
a canard cycle nearby T then it is a hyperbolic repeller canard cycle. Under 
some other conditions on Zq (reversing the directions of Xq and Yo) we can 
derive that such canard cycle is an attractor. 

2.1.3. Analysis around the two— fold singularity. In Equation ([I]), for j3 > 0, 
it is possible to define a First Return Map ip 1 ^ : {\f]5, 3y/j3 + /j,) — > (\/]3, 3\f]3), 
associated to Z\^^, given by 

where £>y M/3 (x) is the first return to X of the orbit— arc of Y^p that passes 
through p = (x,0) and gx x (%) is the first return to £ of the orbit— arc of 
X\ that passes through p = (x , 0) . 

Lemma 9. If (3 > 0, A = y/]3 and fi = in ([I]) rAen ('see Figure^) the First 
Return Map ip^(x) satisfies 

(i) ?p x (x) < x ,\/x £ (Vft 3V 7 ?) and 

(ii) 

Proof. Consider Figure [9l Given a point p S L2, the positive Y"— orbit by p 
reaches L3 at the point g = (qi, #2) and the negative X— orbit by p reaches 
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Lq at the point p = (pi,P2)- The negative Y — orbit by p reaches L3 at the 
point q = (91,92). Since 

_ (pi - 3V^)(pi - VP) 3 (pi(17U + 99V?) - V?(6256 + 321/3)) 
92 " g2 " 384^ 

and y/JJ < pi < 3\ff3 we conclude that 92 — 92 > and item (i) is proved. 
Item (ii) follows from Section T2.1.21 □ 

Note that Lemma [9] implies that Z^-p p does not have closed orbits in the 
interior of the closed curve of Z passing through the visible S— fold point 
of Yob- Moreover, when p < (see Figure [T0|) . Lemma [9] guarantees that 
ipt? has a unique fixed point x where x < 3yf\f3) + p. And, in this case, 
\(il>x)'(%))\ ^ 1> i- e -' ^ is a hyperbolic fixed point for tp 1 ^ that corresponds to 
a hyperbolic canard cycle of Z\p % p. When p > (see Figure [TO]) , tp^(x) < x 
for all x G (\fj3, 3\/]3 + p) and closed orbits of Z\ : p^ do not arise. 





3^/P + p3^ 3y/p y/p 

p<0 p = 

FIGURE 10. Graph of the First Return Map ip{ 



3^3,/p + p 
p>0 



Given Z = (X, Y), we describe some properties of both X = X\ and 
either Y = Yp tfi or Y = Yp. 

The parameter A measures how the X— fold point d = (A, 0) of X is trans- 
lated away from the origin. More specifically, if A < then d is translated 
to the left hand side and if A > then d is translated to the right hand side. 

The parameter /3 distinguishes the contact order between a trajectory of 
Y and S. In this way, it occurs one, and only one, of the following situations: 
• Y + : In this case /3 > 0. So Y has two T,— fold points in such a way 
that one of them invisible and the other one visible. These points are 
expressed by a = ap = (—y/~j3,0) and b = bp = (y/Jf, 0). Moreover, 
a third point c = cp^ = (3\f]3 + p, 0) plays an important role at 
the analysis of (pQ). This point is the locus where the orbit— arc 
7y (a) intersects transversally £ for negative time (see Figure [TT1) . 
Using the bump function B the distance between c and b is bigger 
or smaller than the distance between a and b according to the value 
of the parameter p. This fact will be important to change from 
Theorem 1 to Theorems 2 and 3. 
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• Y°: In this case /3 = 0. So Y has a S— cusp point e = (0,0) (see 
Figure CQ) . 

• Y~: In this case j3 < 0. So Y does not have S— fold points. In this 
way, Y.f ^ and Y is transversal to S (see Figure [T2l). 




Figure 11. Case Y + . Figure 12. CaseF - . 



2.2. The Direction Function. The next function will be very useful in 
the sequel. 

On E, consider the point C = (Ci, C2), the vectors X(C) = (D±, D2) and 
Y(C) = (Ei, E2) (as illustrated in Figure [T3|) . Observe that the straight line 
r(C) by q + -X"(<z) and g + generically, meets £ in a point p(C). We 

define the C r — map 

p : £ — ► £ 
z 1 — ► p(z). 

We choose local coordinates such that £ is the x— axis; so C = (Ci,0) 
and p(C) 6Mx {0} can be identified with points in R. According with this 
identification, the direction function on £ is defined by 

H : R — ► E 

z 1 — >• p(z) — z. 



r(C7); 




c + y(C7) 



/C + X(C) 



FIGURE 13. Direction function. 



Remark 4. We obtain that H is a C — map. When C G E e U E s i/ie 
following holds: 

• if H(C) < i/ien i/ie orientation of in a small neighborhood of 
C is from B to A; 

• i/ ff(C) = tfien C £ TP; 

• if H(C) > then the orientation of in a small neighborhood of 
C is from A to B. 
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Simple calculations show that p{C\) = E2<yDl +C ^) — +Cl - and conse- 
quently, 

(z\ ntr \ - E<2Dl ~ D2El 

(8) H(Cl) ~ E 2 -D 2 ■ 

Remark 5. If X.f(p) = and Y.f(p) ^ then, in a neighborhood V p of 
p in S, holds H{V p )Di > 0, where X{p) = (D 1 ,D 2 ). In fact, X.f(p) = 
and Y.f(p) ^ are equivalent to say that D 2 = and Eo ^ in (j8j) . So, 
lim H(p\) = D\, where kg ^ and p = (pi,p 2 ). 

(D 2 ,E 2 )^{0,k Q ) 

Considering the previous notation and identifying £ with the x— axis, we 
have that r(C) n £ = when £?2 = ^2- m such a case, H is not defined at 
C. The following property is immediate. 

Proposition 10. // n\ is the number of pseudo equilibria and n 2 is the 
number of virtual pseudo equilibria then n\ + n 2 = v\ + v 2 where v\ is 
the number of zeros of H and v 2 is the number of points q of T, such that 
r(q) n £ = 0. 

Proof. Straightforward according to Remark HJ Equation ([8]) and Definition 

H 

□ 

Remark 6. Given Z\^^, we list some properties of the function H . Ac- 
cording to ([8j) we have that the expression of H is 

Hi(x,X,P,n) 



H 2 (x,X,/3,ix) 



dB 

where Hi(x, A, f3, fi) = —x 2 — x + A + j3 — (3, /i) and H 2 (x, A, /?, /i) = 

ox 

dB 

-x 2 + x-\ + /3 — —(x,(3,fi). So, 
ox 

(i) When x = A we get Hi(\, A, /3, fi) = H 2 {\, A, /3, /i). 

dB 

(ii) For the parameter values satisfying (3 = A 2 + (A, /3, n) > we 

have Hi (A, A, /3,/jl) = H 2 (A, A, p, fi) = 0. 

(iii) Since Hi(0, 0, 0, 0) = (respectively H 2 (0, 0,0,0) = 0) and 
OH dH 

—2(0,0,0,0) = -1 (respectively —2(0,0,0,0) = I), by the Im- 
dx ox 
plicit Function Theorem there is a unique x = ^^(A,/?,^) such that 

Hi(x Hl (X,/3,fi),X,P,^) =0 (respectively H 2 (xh 2 (K P, v), A,/3,/x) = 

0). Therefore, there is only one zero of Hi and only one zero of 

H 2 in a sufficiently small neighborhood of x = 0. These points are 

called pi and r\, respectively, in Figure [14\ The pseudo equilibrium 

pi and the virtual pseudo equilibrium ri are the unique roots of Hi 

and H 2 , respectively, that are relevant to our analysis. In fact, the 

other roots are far from the origin. 
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2.3. Proof of Lemma [8j Now we prove Lemma 

Proof of Lemma\M Here we construct a £— preserving homeomorphism h 
that sends orbits of Z = (X,Y) to orbits of Z = (X,Y), where Z = Z p 
is given by © with p\ = 1 and pi = —1, i = 2,3,4. The other choices 
on parameters pi, i = 1,2,3,4, are treated in a similar way. Let p (respec- 
tively, p) be the fold— cusp singularity of Z (respectively, Z) (see Figure [T5l) . 
Identify p with p, i.e., h(p) = p. Consider a point q S 7 (respectively, 967), 




FIGURE 15. Construction of the homeomorphism. 
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where 7 (respectively, 7) is the orbit— arc of Y (respectively, Y) starting at 
p (respectively, p). Identify 7 with 7 (i.e., h(^) = 7) from a reparametriza- 
tion by arc— length. Let T (respectively, T) be transversal sections to Y 
(respectively, Y) passing through q (resp., q) with small amplitude. Identify 
T with T (i.e., h(T) = T) by arc— length. Let ql G T be a point on the 
left of q. Using the Implicit Function Theorem (abbreviated by IFT), there 
exists a time tl < 0, depending on ql, such that <i>y{<ll-,tl) := p\ G S. Since 
h(T) = T, there exists ql G T such that h(q~) = qj. Using IFT, there exists 
a time tl < 0, depending on ql, such that 4>y(ql,tl) := pi G S. Identify 

the orbit— arc a p {{Y) of Y joining p l to ql with the orbit— arc a p ~(Y) of V 

joining pl to (i.e., h(a p f(Y)) = a p ~(Y)) by arc— length. Fix the notation 

for the orbit— arcs of a given vector field joining two points. Since p (re- 
spectively, p) is a S— fold point of X (respectively, X), using the IFT, there 
exists a time t 2 s > (respectively, l? s > 0), depending on (respectively, p^), 
such that 4>x(pl,t 2 ) '■= p 2 s G S (respectively, ^(p],t^) := G S). Identify 

ai(X) with 5^(X) (i.e., /i(ai(X)) = ct~'(X)) by arc-length. Using the 

Ps pl Ps pi 

IFT, there exists a time > (respectively, > 0), depending on (re- 
spectively, such that (f) Y (Ps,t 3 s ) := q 2 s G T (resp., := §f G f ). 

Identify cr 9 !(y) with a q \{Y) (i.e., /i(a 9 l(y)) = o^ff?)) by arc-length. 

Ps p2 p s ^ p 2 

So, the homeomorphism h sends S to S and sends orbits of Z to orbits 
of Z. 

□ 

3. Proof of Theorem 1 

Proof of Theorem 1. In Case li we assume that Y presents the behavior Y~ 
where f3 < 0. In Cases 2\, 3± and 4i we assume that Y presents the behavior 
Y° where (3 = 0. In these cases canard cycles do not arise (for a proof, see 

0)- 

o Case li. P < 0: The points of X on the left of d belong to S e and 
the points on the right of d belong to E c . See Figure [TBI Since j3 < 0, 
the graph of H is illustrated in H-3 of Figure HH We get that p\ = (— 1 + 
V1 + 4/3 + 4A/2, 0) G S e is a S-repeller and r x = (1 - Vl + 4/^-4A/2, 0) G 
S c . 

o Case 2i. A < 0, Case 3i. A = and Case 4i. A > 0: The configuration 
of the connected components of £ is the same as Case li. Since j3 = 0, the 
graph of H, when A ^ 0, is given by H-3 of Figure [TH When A = (Case 
3i), the graph of H is given by H-2 of Figure PT4"1 and p\ = r\. These cases 
are illustrated in Figure [TTl 

In Cases 5i — 17i we assume that Y presents the behavior Y + where 
13 > 0. 
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Figure 16. Case l x . 




o Case 5i. A < — y 7 ?-' The points of E on the left of d belong to E e , the 
points inside the interval (a, b) belong to E s and the points on (d, a) and on 
the right of b belong to E c . The graph of H is like H-3 of Figure [141 We 
can prove that pi is a X— repeller situated on the left of d and n G (d,a). 
canard cycles do not arise. See Figure fl8l 




Figure 18. Case5i. Figure 19. Case 6i. 



o Case 6i. A = In this case the points on the right of b belong to 

E c , the points on (a = d, b) belong to E s and the points on the left of a = d 
belong to S e . Since /3 = A 2 , H is like H-2 of Fi gure [14] and pi = n . There 
exists a non hyperbolic canard cycle T of kind III passing through a and c. 
See Figure [I9j 

o Case 7i- -y 7 ^ < A < 0, Case 8\. A = and Case 9\. < A < ^//3: 
The configuration of the connected components of E is like Case 5i replacing 
a by d and vice- versa. The graph of H is like H-l of Figure [TH We observe 
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that pi G (d,b) is a X— attractor and n G (a,d). There exists a hyperbolic 
repeller canard cycle T of kind III passing through a and c. See Figure 1201 



b < d < 



d = 



< d < c 




Figure 20. Cases 7i - 9i. 



o Case 10i . A = yffi: In this case the points on the left of a belong to 
X 6 and the points on the right of a belong to X c , except by Q = (6, 0) G E. 
Since /3 = A 2 , H is like H-2 of Figure [HI and p\ =r\. Since fi = and d = b, 
by the construction of the bump function B it is straightforward to show 
that the point Q behaves itself like a weak attractor for Z and there exists a 
non hyperbolic canard cycle of kind III passing through a and c. See Figure 
[3l This case has already been discussed previously in Subsection 12.11 Note 
that in [9] the authors avoid this case. 

o Case Hi. sf$ < A < L\: The meaning of the value L\ will be given 
below in this case. The points of S on the left of a and on (b, d) belong to 
S e . The points on (a, b) and on the right of d belong to S c . The graph of 
H is like H-3 of Figure Q31 We can prove that p\ G (b, d) is a S— repeller 
and r\ is on the right of d. Since the point Q of the previous case is a weak 
attractor, in a neighborhood of d occurs a Like Hopf Bifurcation. Moreover, 
according to Lemma [9l there is a unique canard cycle Ti in a neighborhood 
of d and a unique canard cycle T2 in a neighborhood of c. Observe that 
both are of kind I, T\ is attractor, T2 is repeller and T\ is located within the 
region bounded by T2- See Figure[2TJ Note that, as A increases, Ti becomes 
bigger and T2 becomes smaller. When A assumes the limit value L\, one of 
them collides with the other. 

o Case Yl\. A = L\: The distribution of the connected components of X 
and the behavior of H are the same as Case Hi. Since A = L%, as described 
in the previous case, there exists a non hyperbolic canard cycle T of kind I 
which is an attractor for the trajectories inside it and is a repeller for the 
trajectories outside it. See Figure I2T1 

o Case 13i. Li < A < 2yf$, Case 14i. A = 2y/P, Case 15i. 2^ < A < 
3tJ]3, Case 16\. A = 3y/~@ and Case 17\. A > 3y/]3: The distribution of 
the connected components of S and the behavior of H are the same as Case 
Hi. Canard cycles do not arise. See Figure l22l 

The Bifurcation Diagram is illustrated in Figure l23l □ 
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FIGURE 21. Cases Hi and 12 1 . 




15i >^ 


f 









FIGURE 22. Cases 13i - 17i. 




A = 2^ 
A = 3v^ 



FIGURE 23. Bifurcation Diagram of Theorem 1. 



Remark 7. In Cases 9i and Hi the ST— bifurcations (as described in [9], 
Subsections 11.2 and 12.2) arise. In fact, note that the trajectory passing 
through a, in Case 9i, and c, in Case Hi, can make more and more turns 
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around p\ . This fact characterizes a global bifurcation also reached in other 



cases. 



4. Proof of Theorem 2 



Proof of Theorem 2. In Case I2 we assume that Y presents the behavior 
Y~. In Cases 22, 32 and 42 we assume that Y presents the behavior Y°. In 
Cases 52 — 192 we assume that Y presents the behavior Y + . 

o Case 1 2 . P < 0, Case 2 2 . A < 0, Case 3 2 . A = 0, Case 4 2 . A > 0. Case 
52- A < —yf$. Case 62. A = Case 7 2 . < A < and Case 82. 

A = 0: By the choice of the bump function B, these cases are analogous to 
Cases li, 2i, 3i, 4i, 5i, 61, 7i and 81. 

o Case 92- < A < — ///2: The analysis of this case is done in a 
similar way as the Case 9i. In this case and in Cases 72 and 82 there exists 
a hyperbolic repeller canard cycle T of kind III passing through a and c. 

o Case 10 2 . A = v 7 /? — fi/2: The points of E on the left of a belong to 
E e and the points on (d,b) belong to E s . The points on (a,d) and on the 
right of b belong to E c . The graph of H is like H-3 of Figure Q3J Observe 
that pi £ (d, 6) is a S— attractor and n is on the right of b. In this case the 
arc Jx( a ) of AT passing through o returns to E at the point c. So, in this 
case there arises a non hyperbolic canard cycle T = 7x(a) U 7y (c). By the 
discussion on subsection 12.1.21 we have that T is a repeller and we do not 
have other canard cycles inside T. See Figure! 





FIGURE 24. Cases 10 2 - 12 2 



o Case 11 2 . \f]3 — fi/2 < A < \f]5: The configuration on E and the graph 
of H are the same as Case 10 2 . Since Q^{c) G (a, d) there exists a point Q £ 
(c),Qx (b)) such that 77' (Q) = 1. So there exists a hyperbolic repeller 
canard cycle V, of kind I, passing through Q. See Figure [241 Moreover, by 
Lemma [9] this canard cycle is unique. In Figure [TO] we introduce the point 
x which plays the same role of Q. 

o Case 12 2 . A = y/]3: The points of E on the left of a belong to E e and the 
points on the right of a belong to E c , except by the tangential singularity 
c = d. The graph of H is like H-2 of Figure [TH The repeller canard cycle 
r presented in the previous case is persistent. Recall that this canard cycle 
is born from the bifurcation of Case 10 2 . So, the radius of T does not tend 
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to zero when A tends to \f]i. Moreover, the tangential singularity b = d 
behaves itself like a weak attractor. See Figure l24l 

o Case 13 2 . ^ffi < \ < L x , Case 14 2 . A = L ± , Case 15 2 . L\ < A < 
2y/P + fi/2, Case 16 2 . A = 2y^+/V% Case 17 2 . 2^/P + fi/2 < A < 
Case 18 2 . A = 3\ffi + /i and Case 19 2 . A > 3\ff3 + ^: The analysis of these 
cases is done in a similar way as Cases Hi, 12±, 13i, 14i, 15i, 16i and 17i, 
respectively. 




FIGURE 25. Bifurcation Diagram of Theorems 2 and 3. 



The bifurcation diagram is illustrated in Figure [25j □ 

5. Proof of Theorem 3 

Proof of Theorem 3. In Case I3 we assume that Y presents the behavior 
Y~ . In Cases 23, 33 and 43 we assume that Y presents the behavior Y°. In 
Cases 53 — 193 we assume that Y presents the behavior Y + . 

o Case I 3 . /3 < 0, Case 2 3 . A < 0, Case 3 3 . A = 0, Case 4 3 . A > 0, Case 
5 3 . A < -V?, Case 63. A = -y/fi, Case 7 3 . -y/P < A < 0, Case 83. A = 
and Case 9^. < A < y/]3: By the choice of the bump function B, these 
cases are analogous to Cases li, 2±, 3i, 4i, 5i, 61, 7i, 81 and 9i. 

o Case IO3. A = \ffi: The distribution of the connected components 
of S and the behavior of H are the same as Case 12 2 . This case differs 
from Case 12 2 because, as observed in Subsection 12.1.31 when A = yf]3 and 
fj, > canard cycles of Z do not arise (see Figure [TO]) bifurcating from the 
non hyperbolic canard cycle T of Case 123 below. Moreover, the tangential 
singularity d = b behaves itself like a weak attractor. See Figure [26l There 
exists a hyperbolic repeller canard cycle V of kind III passing through a and 
c. 

o Case II3. < A < v 7 /^ — A i /2.' The points of S on the left of a and 
on (6, d) belong to S e . The points on (a, 6) and on the right of d belong to 
S c . The graph of H is like H-3 of Figure [Til We can prove that p\ G (6, d) 
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Figure 26. Cases 10 3 - 12 3 . 



is a £— repeller and r\ is on the right of d. Since Qy{qx{o)) £ (a, b) there 
exists a point Q G (^(wWliSrM) such that r/(Q) = 1. So there exists 
a hyperbolic attractor canard cycle T, of kind I, passing through Q. See 
Figure [26J By Lemma in this Hopf Bifurcation a unique canard cycle 
arises. Moreover, there exists a hyperbolic repeller canard cycle T of kind 
III passing through a and c. 

o Case 123. A = yffi — fi/2: The configuration on £ and the graph of H 
are the same as Case II3. The attractor canard cycle T presented in the 
previous case is persistent. Recall that this canard cycle is born from the 
bifurcation of Case IO3. So, the radius of T does not tend to zero when 
A tends to yffi + fi/2. Moreover, it appears a non hyperbolic canard cycle 
passing through a and c. See Figure [Ml 

o Case 13 3 . y 7 /? - /V 2 < A < Li, Case 14 3 . A = L\, Case 15 3 . L\ < A < 
2y^-/i/2, Case I63. A = 2y/fJ-n/2, Case 17 3 . 2^-^/2 < A < Sy^-fi, 
Case I83. A = 3\//? — /i aric ^ Case 193. A > 3\/^ — M ; The analysis of these 
cases is done in a similar way as Cases Hi, 12i, 13i, 14i, 15i, I61 and 17i, 
respectively. 

The bifurcation diagram is illustrated in Figure [251 replacing the number 
2 subscript by the number 3. □ 

6. Proof of Theorem A 

Proof of Theorem A. Since in Equation (pQ) we can take \i G (— ;Uo, fJ-o), from 
Theorems 1, 2 and 3 we derive that it bifurcation diagram contains all the 55 
cases described in Theorems 1, 2 and 3. But some of them are S— equivalent 
and the number of distinct topological behaviors is 23. Moreover, each 
topological behavior can be represented respectively by the Cases li, 2i, 3i, 
4i, 5i, 61, 7i, 81, 9i, 10i, Hi, 12i, 13i, 14i, 15i, I61, 17i, 10 2 , 11 2 , 12 2 , 10 3 , 
11 3 and 12 3 . 

The full behavior of the three— parameter family of NSDS's expressed by 
Equation ([T]) is illustrated in Figure [27] where we consider a sphere around 
the point (A, f3, fj,) = (0, 0, 0) with a small ray and so we make a stereo- 
graphic projection defined on the entire sphere, except the south pole. Still 
in relation to this figure, the numbers pictured correspond to the occurrence 
of the cases described in the previous theorems. As expected, the cases 3i 
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and 32 are not represented in this figure because they are, respectively, the 
center and the south pole of the sphere. □ 




A < < 1 >A > 

A = 



FIGURE 27. Bifurcation diagram of the invisible fold— cusp singularity. 



7. Proof of Theorem B 

When we consider Equation ([3]) the function H, given by ([8]), is constant 
and equal to 1 independently of the value of [x. Moreover, distinct values 
of the bump function B (where B ^ B) do not produce any topological 
change in the bifurcation diagram of the singularity. In another words, 
two parameters are enough to describe the full behavior of this singularity. 
Observe that, by Proposition [TO] we have = and it does not have 
virtual pseudo equilibria. 
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Proof of Theorem B. Since X has a unique E— fold point which is visible 
we conclude that canard cycles do not arise. In Case Is we assume that 
Y presents the behavior Y~ . In Cases 2b, 3 b and 4^ we assume that Y 
presents the behavior Y°. In Cases 5b — Hb we assume that Y presents the 
behavior Y + . 

o Case 1b- P < 0: The points of E on the left of d belong to E c and the 
points on the right of d belong to E e . See Figure [28l 




o Case 2b- A < 0, Case 3b- A = and Case 4#. A > 0: The configura- 
tion of the connected components of E is the same as Case 1b- Note that, 
when A < (Case 2b), it appears a tangential singularity P = (A,0) 6 E e 
but Z s is always oriented from the left to the right. These cases are illus- 
trated in Figure [29J 



d < e d = e d > e 




Figure 29. Cases 2 B - A B . 



o Case h B - A < -2y//3, Case 6b- A = -2y/J3 and Case 7b- -2yfi3 < A < 
—\f/3: The points of X on the right of b and inside the interval (d, a) belong 
to S e . The points on (a, b) and on the left of d belong to S c . See Figure l30l 

o Case 8b- A = —\f^: In this case a = d and the configuration of the 
connected components of E is illustrated in Figure [3T1 

o Case 9b- —y/P < A < \f]5: The points of S on the right side of b belong 
to S e and the points inside the interval (a, d) belong to T, s . The points on 
(d, b) and on the left of a belong to S c . See Figure l32l 

o Case 10b- A = ^ffi: In this case d = b and the configuration of the 
connected components of E is illustrated in Figure [32j 
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FIGURE 30. Cases 5 B - 7b- 




Figure 31. Case 8 B - 

o Case 11b- A > \fj3: The points of S on the right of d belong to S e 
and the points inside the interval (a,b) belong to T, s . The points on (b,d) 
and on the left of a belong to S c . See Figure 1321 




Figure 32. Cases 9 B -11 B . 

The bifurcation diagram is illustrated in Figure [33j □ 
8. Concluding Remarks 

The results in Section 12 of [9] were revisited and extended in this paper. 
The bifurcation diagram of a three— parameter family of NSDS's presenting 
a fold— cusp singularity is exhibited. In particular it is shown the existence 
of some new interesting global bifurcations around the standard fold— cusp 
singularity expressed by ([5]). Moreover, the simultaneous occurrence of such 
local and global bifurcations indicates how complex is the behavior of this 
singularity. 
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A = -2v^ 




FIGURE 33. Bifurcation Diagram of Theorem B. 
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